Relaxation oscillations, stability, and cavity feedback in a superradiant Raman laser 
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We experimentally study the relaxation oscillations and amplitude stability properties of an optical 
laser operating deep into the bad-cavity regime using a laser-cooled 87 Rb Raman laser. By combining 
measurements of the laser light field with non-demolition measurements of the atomic populations, 
we infer the response of the the gain medium represented by a collective atomic Bloch vector. The 
results are qualitatively explained with a simple model. Measurements and theory are extended to 
include the effect of intermediate repumping states on the closed-loop stability of the oscillator and 
the role of cavity- feedback on stabilizing or enhancing relaxation oscillations. This experimental 
study of the stability of an optical laser operating deep into the bad-cavity regime will guide future 
development of superradiant lasers with ultranarrow linewidths. 
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Optical lasers operating deep in the bad-cavity or 
superradiant regime, in which the cavity linewidth k 
is much larger than the gain bandwidth 7^, have at- 
tracted recent theoretical[TJ[2] and experimental [3 inter- 
est. The interest has been partially driven by the possi- 
bility of creating spectrally narrow lasers with linewidths 
< 1 millihertz and dramatically reduced sensitivity to 
the vibrations that limit state of the art narrow lasers 
and keep them from operating outside the laboratory 
environment [4] . These lasers may improve measurements 
of time[5], gravity [6], and fundamental constants|7J [8] 
aiding the search for physics beyond the standard model. 
The cold-atom Raman superradiant laser utilized here 
operates deep into the bad-cavity regime (k/j± ~ 10 3 ^> 
1), making it an important physics test-bed for funda- 
mental and practical explorations of bad-cavity optical 
lasers. 

In the interest of fundamental science and in light 
of the potential applications, it is important to under- 
stand the impact of external perturbations on lasers op- 
erating deep into the bad cavity regime. In this Let- 
ter, we present an experimental study of the response to 
external perturbations of the amplitude, atomic inver- 
sion, and atomic polarization of an optical laser operat- 
ing deep into the bad-cavity regime. In contrast, exper- 
iments have extensively studied the amplitude stability 
properties of good-cavity lasers (k « 71) (See Ref. [9] 
and references therein). Previous experimental work in 
the extreme bad cavity [3] and crossover regime [10] fo- 
cused on the phase properties of the light and atomic 
medium. Amplitude oscillations, intensity noise, and 
chaotic instabilities have been observed in gas lasers 
operating near the cross-over regime (ft/7± < 1Q)[TH- 
[14]. Relaxation oscillations of the field have been stud- 
ied deep into the bad cavity regime using masers[T5] 
in which the radiation wavelength is comparable to the 
size of the gain medium, unlike in the present opti- 
cal system. Previous theoretical studies of amplitude 
stability deep in the bad-cavity regime include stud- 
ies of relaxation oscillations [16], chaotic instabilities [17], 



and intensity fluctuations characterized by correlation 
functions US EB]. 

In good-cavity optical lasers, the atomic polarization 
(proportional to J±) can be adiabatically eliminated and 
the relaxation oscillations are associated with the flow of 
energy back and forth between the gain inversion (pro- 
portional to J z ) and the cavity field A, where J z , J±_ 
are components of the collective Bloch vector J describ- 
ing the atomic gain medium. In contrast, in a bad- 
cavity laser, the cavity field can be adiabatically elim- 
inated, and the oscillations are driven by the coupling 
of J± and J z . Here, we will measure and infer not only 
the light field A(t), but also the atomic degrees of free- 
dom J±(t) and J z (t) using non-demolition cavity-aided 
measurements [T9] to give the complete picture of the dy- 
namics of relaxation oscillations in a bad-cavity laser. 

We will also consider the effects on the laser's ampli- 
tude stability of non-ideal repumping through multiple 
intermediate states. Intermediate repumping states were 
not included in previous simple theoretical models [T], but 
are present in most actual realizations. In addition, we 
demonstrate that the cavity frequency tuning in response 
to the distribution of atomic population among various 
ground states can be used to suppress or enhance relax- 
ation oscillations. As evidence, we show stabilization of 
J Z1 Jl, and A similar to observations of the suppression 
of relaxation oscillations in good-cavity lasers [20]. The 
cavity frequency tuning mechanism is related to other 
applications of cavity feedback including the creation of 
nonlinearities for generating spin-squeezed atomic ensem- 
bles [21], cavity cooling and amplification in atomic and 
opto-mechanical systems [22] , and the control of instabil- 
ities in gravitational wave detectors [23] . 

Our experimental system consists of a quasi-steady- 
state Raman laser described in Fig. p] and in Ref. [3] . 
The laser uses N = 1 x 10 6 to 2 x 10 6 ^Rb atoms as the 
gain medium. The atoms are trapped and laser cooled 
into the Doppler-insensitive Lamb-Dicke regime (40 /iK) 
in a 1-D optical lattice at 823 nm formed by a standing 
wave in a moderate finesse F ~ 700 optical cavity with 
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FIG. 1. (color online) (a), (b) Physical setup and energy level 
diagram. The trapping light (orange) and Raman dressing 
laser (red, power oc Q? d ) are injected along the cavity axis. 
The repumping light (purple, green) is applied perpendicular 
to the cavity axis. The emitted optical laser light (blue) is 
nearly resonant with the cavity mode (dashed lines) detuned 
from cjcav by 5. The repumping is accomplished through 
a pair of two-photon transitions through intermediate opti- 
cally excited states |II) and | III) with incoherent decay rates 
r. We individually control the two-photon rates W and T3 
with the repumping laser powers oc n? ;2 - |3) represents other 
metastable ground states besides the laser levels, (c) Example 
emitted laser photon flux |A(t)| 2 versus time showing spiking 
and relaxation oscillations at turn-on. 



a cavity power decay rate k)2i\ = 11 MHz. The atom- 
cavity coupling to the relevant TEMoo cavity mode on 
resonance with the emitted light is characterized by the 
single particle cooperativity parameter of cavity-QED 
C = 8 x 1(T 3 < 1. 

Fig. [IJd shows a simplified energy level diagram 
of a three level Raman laser system. The las- 
ing transition is a spontaneous optical Raman tran- 
sition from |t) \5 2 S 1/2 F = 2,ra F = 0) to ||) 
|5 2 5'i/2^ = l,rriF = 0) induced by a dressing laser de- 
tuned from the |t) to |I) = |5 2 Pi/ 2 F' = 2) transition by 
A/27T = +1.1 GHz. The incoherent repumping process 
takes place through a third metastable ground state |3) 
which represents the sum of the other hyperfine ground 
states in 87 Rb. The full energy level diagram with details 
of the dressing and repumping lasers is provided in Ref. 

We control the dynamics of the system primarily 
through three two-photon population transfer rates: 7, 
W, and T3 . We control 7 (typical value 60 s _1 ) using 
the intensity of the 795 nm dressing laser. We control the 
repumping rates W and using two 780 nm repumping 
lasers tuned near resonance with the \5 2 S 1/2 F = 1,2) 
\5 2 Ps/2F' = 2) transitions. The repumping intensities 
are independently controlled allowing us to set the pro- 
portionality factor r = Ts /W . In our experiments, 



W ranges from 10 3 s _1 to 10 5 s _1 , and r ranges from 
0.01 to 2. The repumping dominates all other homoge- 
nous broadening of the |t) to \\) transition such that 
7_l & W/2. The inhomogenous broadening of the tran- 



sition is 7^ 



10 3 



To summarize, the relevant 



rates to characterize our system satisfy the hierarchy 
k > 7± ~ W/2 ~ NCj > j in > 7. Coupling to other 
transverse and longitudinal cavity modes is negligible. 

The frequency of the superradiantly emitted light uo 1 is 
set by the frequency of the dressing laser and the hyper- 
fine splitting cjhf/27t = 6.834 GHz. The detuning of light 
and cavity resonance frequency is 5 = 2(co> cav — cj 7 )/ac, 
normalized to the cavity half linewidth. The single par- 
ticle scattering rate from the dressing laser into the cavity 
modeisr c (£) = C'7/(l + £ 2 ). 

The cavity frequency is dispersively tuned by the 
atomic ensemble cj cav = LUbcav J r^Zt a k^'k, where 

^bcav IS 

the bare cavity frequency and a& is the cavity frequency 
shift for a single atom in the kth ground Zeeman state 
resulting from dispersive phase shifts of the intracavity 
light field. 

Since ^ otj f° r & 7^ j\ the cavity frequency can 
provide a measurement of the populations Nk- We can 
suddenly switch off the repumping and dressing lasers 
to effectively freeze the atomic populations [25], We 
then combine repeated non-demolition cavity frequency 
measurements p~9j [26j [27] and NMR-like rotations [28] to 
determine J z {t) = {\ £f (|t;> (til - \U) <W)> and 5{t) 
[24] . The amplitude of the light field emitted from the 
cavity A(t) just prior to freezing the system, along with 
the measurement of the cavity frequency detuning 5 pro- 



J- 



using the relation 

N 



vides an inferred value of J± - 

A(t) = y/T c (5(t))J ± (t), where J_ = £f I!,) (t 

We observe characteristic laser spiking and relaxation 
oscillation behavior in |A(t)| 2 as the laser turns on and 
settles to steady state (Fig. [l^). To systematically study 
small amplitude deviations about the steady-state val- 
ues of J 2 , s , A s , and Jj_ jS , we apply a swept sine tech- 
inque, similar to Ref. [29]. We apply a simultaneous 
small amplitude modulation of the repumping rates as 
W(t) = W{1 + ecoscjt)) and T 3 (t) = rW(t). The mod- 
ulation frequency uj is scanned over frequencies of order 
W, such that 7 < uj <C n. We then measure and infer 
the quantities A(t), J±(i), and J z (t) as described earlier. 

To calculate the transfer function of the applied mod- 
ulation, the measured light field amplitude A(t) exit- 
ing the cavity as a function of time is fit to A(t) = 
A s (l + a(uj) cos(ujt + (j) a (uj))). The normalized fractional 
amplitude response transfer function is Ta(cj) = a(uo)/e 
and the phase response transfer function is = (j) a (u). 
We also define the modulation frequency that maximizes 
Ta(cj) as the resonance frequency u res . We present the 
results in Figs. [2] [3] and [4] which we discuss in detail 
after developing a theoretical model to guide the inter- 
pet at ion. 
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FIG. 2. (color online) (a) Parametric plots of the Bloch vector 
components J z (£) and J± (£) over a single cycle of modulation 
of the repumping rate W near Wpk for modulation frequen- 
cies below, near, and above resonance or oj/2tv = 0.11,2.2,8 
kHz, from left to right. The black points are the mea- 
sured small-signal deviations about the measured steady-state 
Bloch vector (red arrow). The blue curve is the predicted 
deviation from steady-state given the experimental parame- 
ters N = 1.3 x 10 6 , r = 0.71, S s = 1, W = OANCj, and 
NCj = 125 x 10 3 s _1 . The modulation depth e for data be- 
low and above resonance was doubled to make the response 
more visible, (b) The amplitude (upper) and phase (lower) 
response transfer functions Ta(cj), T^uj) of the light field for 
three values of the repumping rate W . The points are mea- 
sured data, and the lines are zero free parameter predictions 
of the response, (inset) cj re s versus W (points) and a fit to 
LUo (line) showing the expected frequency dependence of the 
relaxation oscillations on repumping rate. 



The theoretical model presented here consists of equa- 
tions for the three ground state levels shown in Fig. [T}}, 
which qualitiatively describes the eight level 8T Rb sys- 
tem. The data is also quantitatively compared to an 
eight-level model described in Ref. [30]. Our theoreti- 
cal model uses the semi-classical optical Bloch equations 
for an effective two level system where n ^> W, 7 so that 
the cavity field is effectively slaved to the atomic po- 
larization and can be adiabatically eliminated from the 
system of equations. Additionally, we have adiabati- 
cally eliminated the populations in the optically excited 
states 1 1) , I II) , I III), arriving at the steady state solutions 
for the inversion J zs and collective atomic coherence 
^±,s[B[3l]. The steady state amplitue A s is maximized 



at W = Wpk = ^NT C (S S ) where 5 S is the steady state 
cavity detuning. 

To predict relaxation oscillations and damping, we 
do a straightforward expansion about the steady state 
values J zs and J± jS as J z (t) ~ J z ,s0- + 3z(t)) and 
J±(t) ~ J±,s0-+j±(t)) and ignore terms that are second 
order in small quantities j±,j z and repumping modula- 
tion amplitude ep~6j[32]. This results in a system of cou- 
pled differential equations for j z ,j± which can be recast 
as two uncoupled higher order differential equations 



P'j'z,± +L,± + 2lo'jz,± + Wo3z,± = D Zj± (uj)ee lu;t (1) 

which, ignoring the third order term, is similar to the 
equation for a damped harmonic oscillator described by 
a damping rate 70, a natural frequency c^o, and a com- 
plex, frequency dependent drive D Zi ±(lj). To discuss the 
stability of the system, we also define the damping ratio 
C = 2^- At critical damping £ = 1, no relaxation oscil- 
lations are present, as the system exponentially relaxes 
to steady state. 

We will consider relevant limiting cases of these equa- 
tions to concisely guide the interpretation of the data. 
The full equations and expressions for 70, cjo, P and 
D Zj ±(lj) are given in Ref. [30] . 

First we consider the limiting case where r — >> 00 and 
S = for which J z and J± respond like a damped 
harmonic oscillator described by 70 = W/2 and ujo = 

However, the drive on 



y/W(NC-f - W) as p = 0. 
the two oscillators is different, leading to different over- 
all responses as D±(uj) = ^(NCj — 2W — iuj) and 
D z (uj) = (NCj — W)(W+ioj). The magnitude and phase 
of the drives change with the modulation frequency and 
repumping rate, even as the modulation depth e remains 
constant. 

We show this driven oscillator response in Fig. [2^i with 
the measured and predicted parametric plot of J z and J±_ 
at three different applied modulation frequencies, with 
repumping near W = W p ^. Although the characteris- 
tic frequencies and rates of the atomic oscillator do not 
change, the differing drives lead to a change in the phase 
relationship between the response of the two quadratures. 
This data demonstrates how the external perturbation 
results in different effective drives for j± and j z , leading 
to different behaviors with modulation frequency. 

In Fig. [2)3, we focus on the light field's transfer func- 
tions Ta(oj) and T^fuj). Data for three different average 
repumping rates W are shown. The data qualitatively 
agrees with Eqn. [TJ and quantitatively agrees with the 
displayed theory calculated for the full model including 
all 87 Rb levels. 

In most experimental systems for a superradiant light 
source, the repumping ratio r is finite. In this case, 
the damping rate and resonance frequencies are 70 = 

r( T(i 2 + 7);l^ 7) and Wo = J^WiNC-y-W). We can 
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FIG. 3. (color online) Effects of finite ratio of repumping rate 
r (a) Comparison at two values of r of Ta(cj) versus modula- 
tion frequency. The points are measured data in good agree- 
ment with the zero- free parameter fit (lines), (b) Plot of the 
resonance frequency cj r es (black) and the peak value Ta (cj re s) 
(red) versus r. 



understand the effect of the nonzero /3 — W (l+ r ^ on the 
response by considering the Fourier transform of Eqn. 
[I] The effective damping rate is the coefficient on the 
sum of the imaginary terms on the left hand side of Eqn. 
[I] 7o = 70 — uj 2 /(W(1 + r)). The frequency dependent 
term associated with f3 lowers the effective damping rate. 
However, as long as 7q > at uo = cjo, the system will 
remain stable. The form of j f also indicates reduced 
damping for r < 1 as shown in Fig. [3^,. We observe an 
increase in the peak relaxation oscillation amplitude as 
r —> (red). From the form of ujo, we also expect to 
see the resonance frequency sharply decrease as r — > 0, 
which we observe in Fig. [SJd (black). 

To understand the dynamic tuning of the cavity reso- 
nance frequency cj cav in response to changes in the atomic 
populations, we next consider the case r — » oo and S ^ 0. 
We also assume the cavity frequency is tuned by the 
atomic inversion J z as a — a± = —a^ > 0. The dy- 
namic cavity tuning then modifies the damping rate as 

7o = \ (l + 2a£ s - c^))' The tunin & of ^cav 

can act as either positive or negative feedback. When 
5 S > the additional term enhances the damping, mak- 
ing the oscillator more stable. If 5 S < 0, the damping can 
become negative, leading to an unstable oscillator. 

We observe this qualitative behavior in the data, which 
agrees well with a full model [30 . As S s is varied from < 
to > 0, the cavity feedback first creates positive feed- 
back that enhances and then negative feedback that de- 
creases the measured amplitude response function Ta(oj) 
in Fig. [4^i. We also measure the effect of the damp- 
ing on J z (t) and J±(t) in Fig. [4J}, where we can see that 
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FIG. 4. (color online) Evidence of negative and positive cavity 
feedback, (a) Amplitude transfer functions of the emitted 
electric field for three detunings from cavity resonance. The 
data points are the average of 4 experimental trials. The lines 
are fitted transfer functions with N as a free parameter, (b) 
Cavity damping of collective atomic degrees of freedom. The 
response going from S s = 0.2 (left) where the 70 is small to 
S s = 0.9 (right) where the system is expected to be critically 
damped. The red lines are sinusoidal fits to the data (circles). 
The dashed lines highlight the damping in both J± and J z . 

both atomic degrees of freedom are stabilized by negative 
cavity- feedback. 

We have studied the dynamics of the polarization, in- 
version, and field of an optical laser operating deep in the 
bad-cavity regime. We have shown that dispersive cav- 
ity frequency tuning can suppress or enhance relaxation 
oscillations. Having experimentally validated our model 
for optical lasers in the extreme bad-cavity regime, future 
work can now extend the formalism to realistic models 
of proposed ultrastable lasers using ultranarrow atomic 
transitions in atoms such as Sr and YbpQ. Further stud- 
ies of the nonlinear dynamics of the extreme bad-cavity 
laser system will include investigations of chaos [17] and 
squeezed light generation [33], 
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FIG. 5. (a) The energy level diagram for the Dl Raman tran- 
sition used for lasing. The linearly polarized Raman dressing 
laser is shown in red and the superradiantly emitted light in 
blue. The cavity mode resonance frequency uj ca v is denoted 
with a blue dashed line. With the quantization axis defined 
by a 2.7 G magnetic field oriented along the cavity axis, the 
linearly polarized light is a linear combination of cr+ and <r_ 
polarizations. (b) The energy level diagram for the D2 re- 
pumping beams F2 (green) and Fl (purple). The dark state 
with respect to the repumping lasers is labeled with a gray 
circle and corresponds to \\). 



Cavity details, full energy level diagram, dressing, 
and repumping laser scheme. 

The optical cavity has a mode waist of 71 /im, mirror 
separation of 1.9 cm, and a finesse ofF = 700. The atoms 
are trapped by a 1-D intracavity optical lattice at 823 
nm and laser cooled to approximately 40 jaK. The sub- 
wavelength localization of the atoms along the cavity-axis 
ensures that the atoms are in the Lamb-Dicke regime 
along this direction. However, the atoms are not in the 
Lamb-Dicke regime with respect to motion transverse to 
the cavity axis. 

We first note that the energy levels, dressing laser, 
repumping lasers, and relative frequency tunings are 
the same as the primary configuration presented in 
Ref. pQ. We include the details here again for clar- 



\b 2 S 1/2 F = l,m F = 0) and 



ity. The ||) 

1 5 2 aSi/ 2 F = 2, vriF = 0) hyperfine ground states of 87 Rb 
form the basis of our Raman laser (energy level dia- 
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gram shown in Fig. |5|. The \\) state is dressed with 
an optical Raman laser at 795 nm that induces an op- 
tical decay to ||), with the cavity tuned to be resonant 
or near-resonant with the emitted light. The intensity 
of the dressing beam and the detuning by A/2tt = +1.1 
GHz from the optically excited intermediate states |I±) = 
1 5 2 Pi/ 2 F' = 2, m' F = ±l) set the single- atom scattering 
rate into all of free space 7. The induced single-atom 
scattering rate into the cavity mode is T c = where 
5 is the detuning of the cavity resonance frequency from 
the emitted light in units of the cavity half-linewidth, 
and C is the single particle cooperativity parameter of 
cavity QED[2 . The branching ratio for decay to ||) are 
included in the definition of C. 

The quantization axis is set by a 2.7 G magnetic field 
along the cavity axis. The Raman dressing laser is in- 
jected non-resonantly along the cavity axis (experimental 
setup shown in Fig. 6a). For this quantization axis, the 
linearly polarized Raman dressing laser is an equal com- 
bination of <r + and cr~ light. Constructive interference 
between the two decay paths from |t) to ||) through the 
two states |I±) leads to enhancement of light emission 
with linear polarization rotated 90° from the polariza- 
tion of the dressing laser. Conversely, the orthogonal 
polarization of emitted light is suppressed by destructive 
interference of the two decay paths. 

Two repumping lasers control the rate out of ||) 
and back into ||). The laser powers are independently 
set using acoustic optic modulators (AOMs) 1 and 2 
shown in Fig. 6a. Both repumpers are 7r-polarized, ap- 
plied perpendicular to the cavity, and separately tuned 
near the |5 2 S 1/2 F = 1, 2) to |5 2 P 3/2 , F' = 2) transi- 
tions. The Fl repumper moves atoms primarily from 
the ground \b 2 Si/ 2 F = l^mp) states to the ground 
|5 2 Si/ 2 F = 2, mp) states, and sets the scattering rate 
VFout of ||). The F2 repumper pushes population 
to \\) as the dipole matrix element for the transition 
|5 2 S 1/2 ,F = 2,m, = 0) |5 2 P3/2,^ = 2,m' f = 0) 
is zero. We quantify the F2 repumping rate by calcu- 
lating the single particle scattering rate for an atom in 
|5 2 5i/2? F = 2, rrif = 1) state as T3 . We define the pump 
ratio r = ^ to quantify to what degree population can 
build up in energy levels outside of ||) and ||), with the 
ideal repumping case being r — > 00. Population will still 
remain in the states |5 2 5'i/2,P = = ±1) even as 

r — » 00, unlike in the ideal three- level model. This is the 
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FIG. 6. Measurement setup and sequence timing diagram, (a) The physical setup for measuring the light field amplitude 
response A(t) and the atomic responses J z (t), J± to a small modulation of the repumping rate W(t). Optical beams are 
colored with arrows; RF and microwave signals are in black, (b) Timing diagram for measurements. Superradiance is started 
at ti, and the repumper W(t) is modulated (with r constant) starting at time £2- At £3, the repumper and dressing lasers 
are shut off, freezing the atomic populations. The microwave 7r-pulse corresponds to a 12 {is pulse of 6.834 GHz microwaves 
resonant only with the ground hyperfine state transition |t) to |^). The pulse completely swaps the populations and N± of 
states |t) and |^). The cavity frequency probe windows Ml and M2 used to determine the dressed cavity frequency oj C av,i and 
w C av,2, and from which we can determine J z {t^). 



reason that the output power correction factor R does 
not asymptote to 1 as r — » 00. Also note that unlike the 
simple three-level model, here, W includes Rayleigh scat- 
tering back to ||), contributing damping of the coherence 
Jj_ without affecting the inversion J z . 

Response function measurement sequence 

In the main text, we measure three quantities to char- 
acterize the response of the superradiant laser to small 
perturbations: the light field amplitude A(t), the mag- 
nitude of the coherence J±(i), and the inversion J z (t). 
As shown in Fig. 6a and b, the perturbation is cre- 
ated by a small fractional modulation of the Fl and 
F2 repumping laser power at frequency uj such that 
W(t) = W(l + e cos ujt) while r(t) remains constant. At 
time £3, we shut off the repumping and dressing lasers in 
less than 100 ns, freezing the atomic populations. The 
amplitude of the light just before shut off A(ts) is deter- 
mined from the IQ-demodulated heterodyne signal shown 
in Fig. 6a. 

We then use methods similar to those in [2] to non- 
destructively probe the cavity-mode to determine the 
cavity resonance frequency uJ ca v,i and hence detuning 
from the emitted light frequency S(ts). We then cal- 
culate the coherence J±(t$) using the relation A(ts) = 
^i(fe)vCV (1 + ^ 2 (ts))- The probe is non-destructive 
in that a small fraction of the atoms are lost or Raman 
scattered to other states during the measurement. 

The inversion J z (ts) is determined by using a mi- 
crowave 7r-pulse to swap the populations between |t) and 
||), and then measuring the cavity resonance frequency 
a second time uJcav,2- Shifts in the cavity frequency due 
to atoms in other states are common mode to both mea- 
surements, such that the difference between the two cav- 
ity frequency measurements is only proportional to the 
inversion uj ca v,i — ^cav,2 = 2(a^ — aqj J z . The frequency 



tuning of the cavity resonance per atom in each state a^^ 
is calculated from the known cavity geometry, atomic 
dipole moments, and bare cavity detuning from reso- 
nance with nearby optical atomic transitions at 795 nm. 
The process is repeated for different stopping times £3 in 
order to sample the modulation period for a set of times 
U. The data J Z (U) versus J±(U) are shown in several 
parametric plots in the main text. 

In the case of the field transfer functions Ta(uj) and 
T^(cj), continuous time traces A(t) are averaged over sev- 
eral trials at the same modulation frequency and phase to 
enhance signal to noise. The average response is then fit 
to A(t) = A s (l+a(uj) cos (uj t a (^)), with transfer func- 
tions calculated as Ta(uj) = a(uu)/e and T^uo) = §a (cj). 
The measurement is then repeated for different uj. 
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